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Abstract. We study nonlinear transport for two coupled one-dimensional quantum wires or carbon nan- 
otubes described by Luttinger liquid theory. Transport properties are shown to crucially depend on the 
contact length Lc- For a special interaction strength, the problem can be solved analytically for arbitrary 
Lc- For point-like contacts and strong interactions, a qualitatively different picture compared to a Fermi 
liquid emerges, characterized by zero-bias anomalies and strong dependence on the applied cross voltage. 
In addition, pronounced Coulomb drag phenomena are important for extended contacts. 
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1 Introduction 

Transport in interacting one-dimensional (ID) quantum 
wires (QW) has attracted ever-increasing attention over 
the past decade. This interest was sparked mainly by the 
discovery of novel ID materials besides standard (semi- 
conductor or organic chain molecule) systems, such as 
] edge states in fractional quantum Hall bars or carbon nan- 
: otubes. Furthermore, ID QWs are predicted to behave as 
a Luttinger liquid (LL) due to electron-electron interac- 
tions Recent nonlinear transport experiments for 
individual nanotubes have indeed demonstrated iinpres- 
I sive agreement with the LL theory of nanotubes 0]. In 
. these experiments, transport was limited either by the 
' contact resistance to the leads, or by a tunnel junction 
I ("topological kink") within the nanotube. In both cases, 
■ the observed power laws in the (nonlinear) conductance 
have allowed for a consistent explanation in terms of LL 
theory. 

Different transport experiments built up of at least 
i two nanotubes can reveal even more dramatic deviations 
< from Fermi liquid transport. The theoretical predictions 
of Ref. [Q for crossed nanotubes (which are coupled in a 
pointlike way) were recently observed experimentally by 
Kim et al. 0. For longer contacts between the nanotubes, 
Coulomb drag [|| is expected to play an important role 
in addition to the crossed nanotube scenario. Coulomb 
drag can be very pronounced in one dimension and leads 
to quite rich physics. In this paper, we study in detail 
two nanotubes arranged parallel to or crossing each other, 
and briefly discuss more complex setups. Notably, such 
experiments are feasible using present-day technology [0, 

The main part of the paper focuses on the schematic 
geometry shown in Figure yj, where transport through two 



clean (ballistic) nanotubes biased by voltages is stud- 
ied. The nanotubes are brought to contact by crossing 
them under an angle Q. By varying this angle, the effec- 
tive contact length Lc can be changed. We shall address 
the crossover from a point-like crossing, where Lc ~ a 
with the lattice spacing a, to an extended coupling, where 
a Lc < L with the tube length L. Besides our pre- 
vious paper [Q, theoretical predictions for transport in 
such a geometry have been given by other authors [^,|^, 
Their results were largely obtained in the linear 
regime, or focus on either very small or very large Lc- In 
this paper, we cover the full crossover from a point-like 
to an extended contact, and explicitly compute nonlinear 
current-voltage relations. For simplicity, we consider the 
same interaction strength parameter g and Fermi velocity 
vp for both QWs. The LL parameter g equals unity for 
a Fermi gas, and becomes smaller for strong repulsive in- 
teractions. Provided one works on an insulating substrate, 
the LL parameter has only a weak logarithmic dependence 
on the tube length L, and the experimentally observed 
value g « 0.25 indicates strong non-Fermi liquid 

behavior. 

For repulsive interactions {g < 1), the most important 
coupling between the two nanotubes is of electrostatic ori- 
gin [Q . It is responsible for the most relevant operator un- 
der the renormalization group (RG), and can cause perfect 
Coulomb drag at low temperatures. In addition, as indi- 
cated by the data of Ref. , electron tunneling between 
the QWs can be important. Provided one has sufficiently 
strong interactions, tunneling into a LL is irrelevant in the 
RG sense. In that case, one can treat it perturbatively, at 
least for not exceedingly large tunneling amplitude. No- 
tice that otherwise our assumption of clean wires breaks 
down in any case, since good mechanical contact of the 
QWs indicates impurity formation within each QW. 
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Fig. 1. Two crossed nanotubes in contact to reservoirs for ap- 
plied voltage Ui and U2, respectively (schematic view). Differ- 
ent crossing angles O imply different effective contact lengths 



The outline of our paper is as follows. In Sec. ^, we 
briefly review the LL concept and discuss how the the- 
ory can incorporate adiabatically coupled voltage reser- 
voirs in terms of Sommerfeld-like radiative boundary con- 
ditions Sec. H presents a detailed discussion of trans- 
port for the setup in Fig. |l| under the assumption of neg- 
ligible tunneling. In particular, for the special interaction 
strength g ~ 1/2, the full temperature-dependent trans- 
port problem is solved for arbitrary contact length Lc- 
The role of adiabaticity of the contacts to the reservoirs 
for the described effects is addressed in Sec. |[ where we 
study weakly contacted reservoirs in the tunneling limit. 
In Sec. 0, effects of inter-wire tunneling are discussed, and 
in Sec.^ some conclusions and possible applications are 
outlined. 

For clarity, we focus on spinless single-channel QWs. 
The modifications arising for spin-1/2 electrons or carbon 
nanotubes are then straightforward. In addition, the rele- 
vant energy scale is supposed to exceed vp/L, so that we 
can effectively put L — s- co. Below we put h = e = ks ~ ^ 
and V = vp/ g = 



2 Model 

The LL concept for ID metals is most transparent in the 
bosonization representation The electron field oper- 
ator is written as superposition of left- and right-moving 
(p = R/L = ±) fermions ijpa{x) for QW a = 1,2. The lat- 
ter are expressed in terms of canonically conjugate bosonic 
fields Oa{x) and (j)a{x), 

ippaix) ^ {2'!ra)^^^^ exp[~ipkFax] (1) 
X exp[-ip{ngy/'^ea{x) - {{t:/ g)^/^(j)aix)] , 

where a is a lattice constant. The Fermi momentum kpa in 
the two QWs can be made different by uniformly shifting 
the chemical potentials of both reservoirs attached to one 
wire. The kpa sue defined such that the chemical poten- 
tials are ±C/i^2/2 as indicated in Figure [|. Using Eq. (|l|), 



the density operator Pa{x) is 



2kpaX 



(2) 



The uncoupled clean QWs correspond to the standard LL 
Hamiltonian, 



H. = 



(3) 



Let us next address how the coupling of each wire to 
the voltage reservoirs can be taken into account. For adi- 
abatic coupling, the electron densities near the end of the 
QW obey the radiative boundary conditions 

^^^p^„(TV2) + ^^PLa(TV2) = ±^, (4) 
where the p = R/ L = ± moving densities in wire a are 



Ppaix) = (47rg) ^^^dx[g0a +P(/>a] 



(5) 



These boundary conditions need to be enforced for the sta- 
tionary expectation values of the densities ppa in the QW 
near the respective contact to the leads. They hold for ar- 
bitrary impurity scattering within each wire and therefore 
also in the presence of coupling between the two wires. 

In the absence of inter- wire tunneling, the charge cur- 
rent 

Ux) = e{g/7ry/^dte^{x) (6) 

is conserved and independent of x. In the presence of tun- 
neling, however, we need to distinguish Ia{x < —Lc/2) 
and laix > Lc/2). Postponing the discussion of inter- wire 
tunneling to Sec. ^, the conductance of each (impurity- 
free) QW is thus Go = e'^/h in the absence of electrostatic 
inter-wire coupling. The latter then implies a reduction of 
the conductance 

Ga = Ic/Ua , (7) 

since the transport-carrying density waves drag each other, 
which makes them "heavier." For low energy scales, the 
electrostatic coupling can be expressed as a local prod- 
uct of densities in the two wires. To justify the locality 
of inter- wire interactions, one can employ the same rea- 
soning as for the locality of the intra- wire interaction 0. 
Neglecting momentum-non-conserving terms, we obtain 



dxC,{x)pi{x)p2{x) 



(8) 



Via / dxC{x)dx0idJ 



lb 



X sm 



dx(^{x) sin 



2kFix 



2kF2X + ^/^TTg92{x) 



where Via = gVi/n and Vib — Vikpikp2/''r^ ■ The function 
((x) specifies the spatial dependence of the inter- wire cou- 
pling, and for practical purposes, we consider ({x) — 1 for 
< Lc/2 and zero otherwise. 



Cix) = [e{x + Lc/2) - 0{x - Lc/2)] 



(9) 
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In addition, tunneling leads to 



Ho = Vo 



(10) 



pp' 



describing electron transfer between the QWs. The bo- 
sonized form of H2 can be found in Ref. Q|. Wc shall 
turn to a discussion of tunneling in Sec. ^ Josephson-type 
couplings fl^ are irrelevant for g < 1 and hence will be 
ignored in the following. 



3 Electrostatically coupled nanotubes 

In this section, we assume that tunneling can be neglected, 
and consider a system described by the Hamiltonian H = 
Hq + Hi under the boundary conditions (Q). Below we 
separately discuss three cases, namely (a) a strictly lo- 
cal contact with Lc 0, (b) a short but finite contact 
length Lc ~ a, and (c) for the special interaction strength 
g = 1/2, we present an analytical solution valid for arbi- 
trarily long contacts. Of particular interest is the differen- 
tial conductance matrix, 



Gaa' = dIa/dUa' 



(11) 



The off-diagonal conductance G12 (or G21), the so-called 
"transconductance" , is the appropriate quantity measur- 
ing Coulomb drag p. For Ui — U2 — 0, the matrix Gaa' 
describes the linear conductances. Generally, the diago- 
nal conductance Gaa and the conductance Ga defined in 
Eq. (m) show qualitatively the same behavior. 



3.1 Strictly local contact 

Let us start with an ideal point-like contact, where from 
Eq. (^, we get for Lc — *■ the result ({x) = LcS{x). In 
that case, tunneling is always irrelevant for g < 1, and the 
only relevant coupling corresponds to the scaling field Vu 
in Eq. (||), provided g < 1/2. For g > 1/2, the effects of 
tunneling (V2) and of the electrostatic coupling (Vi;,) can 
both be treated perturbatively. We shall therefore focus 
on the most interesting strong- interaction region g < 1/2 
in this subsection, and omit the irrelevant perturbation 
Via as well as tunneling (V2). Introducing symmetric and 
antisymmetric fields [Q, 



e±ix) = [0i{x)±02{x)]/V2 , 



(12) 



the Hamiltonian decouples, H = H+ + H^, with 

H± = l j dx[{d^<p±f + {d^e±f] 

± (Lc^^ib/2)cos[v/8^e±(0)] . 
An effective coupling strength is defined as 

Tb - (c<,/a)[aLcV^ib] 1/(1-23) ^ 



(13) 



where Cg is a numerical constant of order unity |14|,|15[. 
The boundary conditions (^) also decouple in the sym- 
metric/antisymmetric {r — ±) channels. Effective right- 
and left-moving [p = R/ L — ±) densities Ppr(x) for these 
channels can be defined in analogy to Eq. (0). The new 
densities again obey the boundary conditions M), but with 
the effective voltages 



C/i 



1,2 



Ur=± = {Ui ± U2)/V2 



(15) 



It is also useful to define the effective current in channel 
r = ±, see Eq. (||). The current in QW a = 1,2 is then 
given by la = (/+ ± I-)/V2. 

Notably, the full nonlinear correlated transport prob- 
lem of crossed LLs therefore completely decouples into 
two effective single-impurity problems r = ± character- 
ized by effective impurity strength iTg, applied voltage 
{Ui ± U2)/V2, and interaction strength 217. This single- 
impurity problem has been studied in detail, e.g. by Kane 
and Fisher ||l^ , and the exact solution for arbitrary in- 
teraction strength has been given in Ref. [ p5[ by combin- 
ing Eq. (|^) and powerful methods from boundary confor- 
mal field theory. This solution can then be immediately 
applied to the crossed LL transport problem. Below we 
discuss the salient features for the special case g = 1/4. 
These features are characteristic for the strong-interaction 
regime g < 1/2. For weak interactions, g > 1/2, all inter- 
wire couplings may be treated in perturbation theory, and 
results can be found, e.g. in Ref. [||. 

The currents through QW a = 1, 2 are 



la 



ieyh)[Ua-{V+±V.)/V2] 



where V± obeys the self-consistency relation 

1 , T B+iiU±-V±/2) 
2tiT 



V± = 27^5 Imif ( - 



(16) 



(17) 



with the digamma function Similar but more com- 
plicated self-consistency equations need to be solved for 
g 7^ 1/4. From Eq. (|lq), one can verify that the conduc- 
tance matrix (^l]) fulfills the bounds 



< Gn/Go < 1 



1/2 < G12/G0 < 1/2 , (18) 



with Go 



^/h. The current- voltage relation ( [16|) or. 



(14) 



equivalently, the nonlinear conductances Gaa' are very d_if- 
ferent from the corresponding results for Fermi liquids 
The correlation effects are most pronounced for T = 
where perfect zero-bias anomalies, a strong dependence 
of Gil on the applied cross voltage U2, and minima in 
Gil for \Ui\ = IC/2I are predicted. Such effects are distinct 
and dramatic signatures of correlations in a LL, and have 
found evidence in recent experiments on crossed multi-wall 
nanotubes [^. Note that for a Fermi liquid, Gn would nei- 
ther depend on J7i nor on U2. Thermal fiuctuations tend 
to smear out these phenomena, see Fig. ^, but they remain 
clearly discernible. 

In the linear transport regime, |J7i,2| ^ T, by Taylor 
expanding Eq. ([l7| ) in the small parameter {Ui ± U2)/T, 
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we find for (7 = 1/4 the voltage-independent linear con- 
ductance, 



Gi/Go 



l + c<f'(i + c) 



c^Tb/2ttT. (19) 



More interesting is the transconductance, which reads to 
leading order in {Ui ± U2)/T: 



Gi2/Go = ^c3tf''"(l/2 + c) 



(20) 



Clearly, for very small applied voltages (either Ui or U2 ap- 
proaching zero), the transconductance vanishes. As shown 
in Sec. 3.2, the vanishing linear transconductance is a gen- 
eral consequence of the point-like nature assumed for the 
contact. For an extended contact, the linear transconduc- 
tance need not vanish. For low temperatures, |C/i 2I ^ 
T < Tb, Eq. gives G12/G0 ~ 2U1U2/TI, while in 
the high-temperature limit, T 3> Tb, we find G12/G0 — 
TrUiU2TB/8T\ 

Finally, for very strong interactions, g < 0.1, bistabil- 
ity effects were reported in Ref. As these phenom- 
ena require spinless electrons and cannot be observed in 
nanotubes, however, we do not further discuss them here. 
Based on the mapping of the transport problem of crossed 
LLs to two decoupled single- impurity problems, it is never- 
theless straightforward to obtain the conductance matrix 
in closed form in the bistability regime. 

3.2 Short contact 

Next we turn to a short but nonlocal coupling, where Lc is 
of the order of a few lattice spacings a. For short contact 
length, the bosonic fields can be expanded in powers of x 
around the center of the coupling region, x = 0. The sub- 
sequent RG analysis shows that the only relevant coupling 



term is still due to the scaling field Vib, and one arrives 



again at the Hamiltonian 



replacement Vi 



lb 



wit 



The only difference is the 



dx(^{x) cos[2(kpi ± kp2)x] 



(21) 



leading to couplings as in Eq. (|lj). For large Lc, 
<C Tg due to the oscillatory integrand. The result- 
ing physics is similar to the point-like case of Sec. B.ll 



but in addition exhibits linear Coulomb drag te,|9[|. No- 
tably, we can exploit the exact solution of Ref. ||l5| for the 
single-impurity problem to solve this problem too. 

In the remainder of this subsection, we again focus 
on the case of g — 1/4 where the algebra becomes quite 
simple. In particular, one has to solve the self-consistency 
equations ( p7| ) with Tb T^. For clarity, we consider the 
special case U2 = 0, and compute the transconductance 
G21 defined in Eq. (11). Notably, even though U2 = 0, the 
current I2 and hence 6*21 can be finite (Coulomb drag). 
The transconductance follows from 



G2i(C/i,T)/Go = 



1 



V2dU, 



■{V--V+) 



(22) 



where the V± are determined as the U2 — solutions of 
Eq. (|l^ with Tb->T- 



B ■ 



The linear transconductance at T > can now be finite 
due to the nonlocality of the contact [p|. We find 



G2i/G( 



1 1 

= 25±r 



c±r(l + e±) 
c±iP"(i + c±) 



(23) 



where c± = Tg /2itT. The temperature dependence of the 
(7=1/4 linear transconductance ( p3| ) is shown in Fig. 
We observe that for a long contact, Tg ^ T^, the per- 
fectly quantized transconductance e^/2/i is approached at 
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low but finite temperature. This is the absolute drag effect 
reported by Nazarov and Averin for a long extended 
contact. Equation ( |2^ ) then describes how the absolute 
drag effect is thermally destroyed at high temperatures. 

However, as T — > 0, the linear transconductance van- 
ishes throughout the strong interaction regime g < 1/2. 
(This is an artefact of the locality assumption and will 
not be true for extended contacts, see below.) In the zero- 
temperature limit, one can obtain the nonlinear g ~ 1/4 
transconductance for |C/i| ^ Tg and 1/2 = from Eqs. (g3) 
and (llT 



G2i/Go-(C/i/4)2[(r^ 



(24) 

Finally, let us address the role of a difference in the 
Fermi momenta, Skp = kpi — kF2^ which could experimen- 
tally be tuned by varying the mean chemical potential in 
one QW relative to the other. In the absence of electron 
tunneling between the QWs, the nonlinear conductance 
matrix Gaa'{Ui, U2) for a short contact does only depend 
on Skp via the following from Eq. (pl|). In fact, for a 
strictly local contact {T^ = Tb), there is no dependence 
on Skp at all unless there is tunneling. As discussed in 
Sec. 1^, this effect can be used to experimentally disentan- 
gle the effects of tunneling and electrostatic coupling. 



3.3 Arbitrary contact length 

Next we study an arbitrarily long contact, where Lc can 
even approach the system length L. Here tunneling can 
be a relevant perturbation under the RG transformation, 
provided the interaction is not too strong. Following the 
analysis of Refs. for g > g' = V2 - 1 ~ 0.414, 

tunneling is relevant. In fact, for g > g" = 1/^/3 ~ 0.577, 
tunneling (V2) is even more relevant than the electrostatic 
coupling Vib. In this regime, one should then first treat the 
tunneling. Below we focus on the strong-interaction case, 
g < g" . For g' < g < g" , the main effect of tunneling is 
a renormalization of the electrostatic coupling ||l^ . In the 
following we assume that Vib contains this renormaliza- 
tion and then omit the scaling field V2. Since for a short 
contact, tunneling is always irrelevant, such a reasoning 
can be applied for arbitrary Lc in the regime g < g" , in 
particular for the case g = 1/2 investigated below. 

Taking into account the Via operator in the bulk (i.e., 
for Lc = L) causes a renormalization of the interaction 
constants, thereby splitting them [ p^ , 

g^9± = -lf=f= ■ (25) 

For Lc <C i, however, the scaling field V-\ „ is irrelevant. 
Below we shall neglect the weak splitting (p5|). Following 
our analysis, this could only create a problem for very long 
contacts, Lc « L. The case g = 1/2 then permits a full 
solution of this transport problem for arbitrary Lc- The 
resulting effective Hamiltonian is 



H = Ho + Vib / dxC{x) sin 2kFix + \p2ji:di{x) 



X sm 



-2k 



F2X 



^e2i-x) 



where the spatial coordinate along QW a = 2 has been 
changed to —x. As we shall see later, this mirroring of the 
axis is crucial to ensure correct anticommutation relations 
between new fermionic fields. Expanding the product of 
sin terms, cos functions of the sum and the difference of 
the fields 601 emerge. In Ref. ||^, only the difference term 
was kept, since the other term is highly oscillatory for long 
contacts and does not contribute for Lc S> a. However, a 
proper description of the crossover from a short to a long 
contact requires to consider the full coupling (p6|). Below 
we take Eq. (0) for C{x). 



3.3.1 Refermionization 

Remarkably, the transport problem posed by Eqs. (|2£ 
and (^) can be solved exactly for arbitrary Lc by refermion- 
ization. Switching to the chiral fields, 

Vr,l = V^i^Pc ± Oc) , (27) 
we obtain with Skp = kpi — kF2 and 2kp = kpi -\- kp2 



„ — 1 O 



Vu f dxCix)! Y^e 
I P=± 



-ip{vl(x)+ip''R{-x))/V2 I 



We then employ the slightly modified refermionization 
transformation of Luther and Emery p8| , 



_ Vl,2 i2kFix+i(v],(xnv>l(-x))/V2 



(28) 



1^/34(2;) = Jlh^- ^±^2kF2XT^{<pl{-x)T<plix))/^/2 

v27ra 



Special Majorana fermions r]i ensure the correct anticom- 
mutation relations between new operators, see Ref. 
However, all chemical potentials of the new particles can 
be incorporated into the refermionized Hamiltonian (A = 
27TaVu) 

H = ~ j dx[ J2 + 2kpi)^^ix) 

- E 'I']{x){id^T'ikp2W^{x)] 

+ i\ j dx(:{x)^^l{x)^^{x) + ^^[x)^l{x) (29) 
,{x)^l{x)-^l{x)^^(x)] . 



(26) 



Then they disappear in the definition (pq). This prob- 
lem permits an exact solution, since we are left with a 
quadratic Hamiltonian expressed in terms of the ^i{x) 
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operators. We obtain this solution via the equations of 
motion 

{dt±d,)wl^{x,t) = ±XC{x)wl,{x,t) (30) 

(dt ± d^pl^ix.t) = TXCix)^l2{x,t) . (31) 

They describe free chiral fermions with hnear dispersion 
outside of the contact area. The contact region, |a;| < 
Lc/2, then acts as a scatterer. Denoting by aj(fc) (&|(fc)) 
the momentum-space creation operator for fermions mov- 
ing towards (away from) the scatterer, their respective 
particle densities can be related by a transmission matrix 
Aj(fc), 



{■.bUk)bJk) 



i=l 



(32) 



The matrix Dij ca n be f ound in closed form from Eqs. 
and (P), see Sec. IsXl 



To solve the full transport problem, we also need to 
re-express the boundary conditions (^ in the new fermion 
basis. Defining pi{x) — (: tf//(a;)tf'i(x) :) as density of the 
new fermions, for g — 1/2, we get the relation to the 
previous densities ppo^ {x) (where p = ± = R/L): 

Ppi{x) = (1 + 2p)[pi{x) + p2ix)] + (1 - 2p)[p3ix) - pi{x)] 
Pp2{-x) = {2p - l)[pi{x) - p2{x)] - (1 + 2p)[p3{x) + pi{x)] 

The normal ordering should be performed with respect 
to the ground state, given by a Fermi sea filled up to 
k — 2kpi and k = ^2kp2 for channels (1,2) and (3,4), re- 
spectively. Plugging these last relations into Eq. (Q) gives 
the boundary conditions: 



3.3.2 Transmission matrix 

The transmission matrix contains only one independent 
element. Obviously, D13 = D31, Dn = D33, D24 = £'42, 
and D22 = D44 because of the system symmetry. All other 
matrix elements vanish since the channels (1,3) and (2,4) 
fully decouple. Furthermore, since the Hamiltonian (p9) 
conserves the net particle numbers in channels (1,3) and 
(2,4), one has the additional relations Di3(fc) — 1 — -Dii(fc) 
and D24{k) = 1 — D22{k). The last simplification stems 
from the symmetry of the equations of motion and reads 
D{k)=D22{k)^Dn{k). 

Let us now explain how to find D{k). Since a right- 
mover in channel 1 is scattered to a left-mover in chan- 
nel 3 within the contact, we can regard them as a single 
species which is backscattered by the contact. Then the 
problem reduces to the determination of the penetration 
coefficient of ID fermions through a rectangular barrier. 
The corresponding result can be found, e.g. in Ref. M, 



D{k) = 



4kW 



A4|sin[L,x]P+4P|x| 



(37) 



where ^ k^ — \^. Next we insert Eqs. ( ^5| ) and ( ^ ) into 
the boundary conditions ( ^3|) and compute the current 
from Eq. iM). 



3Pi 



-^3p2 +P3 - Pi = Ui/A-K , 
p+ - 3p+ - 3p+ = -(72/4^ , 
p+ + 3p+ - 3p+ = -t/i/4^ , 
'riP2 - P3 ~ Pi ^ C/2/47r , 
where pf = pi{±L/2). For the current, one gets 
h = 2[pi{x) + p2{x) - p3{x) + pi{x)] . 



(33) 



pX- 

-3pr 



3.3.3 Linear transconductance 

The linear transconductance G12 can be found in closed 
form, where we focus on T = and Skp = again. No- 
tably, this quantity does not vanish in general as it would 
for a point-like coupling. Linearizing Eq. ( ^6|) in the chem- 
ical potentials Ui, we need to solve the now linear system 
of equations ( p3[) for the pi, yielding 



2G12/G0 — 



1 - D{2kF) 
l-D{2kF)/2 



(38) 



(34) 



Here x is arbitrary due to the continuity equation. The 
incoming free fermions must obey the Fermi distribution. 



{■.a\{k)aAk) 



np{k) 



--jy-^-jv-j •/ - npik- pj) 
with npiE) = l/[exp (£■/£) + 1]. Therefore we obtain 



For uncoupled nanotubes, A ^ 0, we find D{2kF) = 1 and 
thus G12 = 0. On the other hand, for strongly coupled 
tubes, A — > cx), we get D{2kp) = and hence recover the 
absolute Coulomb drag, G12 = /2h, in this limit. For 
small contact lengths, it is possible to perform a Taylor 
expansion of D{k). The resulting linear transconductance 
behaves according to 



dk 



p1= I ^{■■a]ik)a^ik) ■.)^p, 



(35) 



2Gi2/Go = l{XL,f(X/kFf 



where — sign should be taken for channels 1 and 2, and 
-I- for 3 and 4. The effective chemical potentials pj have 
to be computed self-consistently, see below. The outgoing 
densities are then given by 



The local approach of Section 3.2 for g = 1/2 gives 
2G12/G0 = i(2^1/ifcL,)2 (kpL.f , 



(39) 



(40) 



pf 



1 1 



npik ~ Pi) ~ npik)] , (36) 



where Vib denotes the coupling strength of the related 
local problem, see Eq. (pl|). These two results suggest the 
following relation between the coupling constants 



where ± apply to the channels (1,2) and (3,4), respec- 
tively. 



1 



47ra V 2 



3\^/V A\^ 1 



kp 



(41) 
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Fig. 4. Nonlinear conductance Gi (in units of Go = /h) as 
a function of Ui for g — 1/2, T = 0, A = 4, 172 = 25, and various 
contact lengths. The solid curve is obtained analytically for the 
shortest contact using the point-like coupling (|l3). 



which establishes a connection between the local approxi- 
mation and the exact solution for the system with a short 
contact. The corresponding amplitudes Vj^ are then given 
by a Taylor expansion of Eqs. ( pi] ) to second order in Lc/a, 
implying that = Vu and = Vib{l - 2{kFL^f /i). 

Moreover, employing the Sommerfeld expansion, the 
finite-temperature corrections 



to Eq. 
Z\Gi2(T)/Go = 



Z\Gi2(r) = Gi2(T)-Gi2(r = 0) 
) are of the form 

^2y2 Q2j^Q. ^ 2kF) 1 



6 



{2-D{2kF)Y 



In the general nonlinear case, the current has to be 
evaluated numerically, see Fig. § AU curves show a sin- 
gularity at Ui — > 0, independent of Lc- This singularity 
implies non- vanishing current for zero applied voltage and 
reflects Coulomb drag. 



3.3.4 Absolute Coulomb drag 

The absolute Coulomb drag reported in Refs. |8[^ can 
be recovered by our g — 1/2 solution. For a very long 
contact [Lc — L), H can be diagonalized by means of a 
Bogoliubov transformation. Consider, e.g. channels j = I 
and 3, with 

dk 



Akx 



27r 



Cj{k) . 



Allowing for 5kp ^ Q and switching to 

ai{k) — ci(fc) cos(5fc + ic3{k) sinSk , 
asik) — — ci(fc) sinSk + ic3{k) cosSk 



where 2dk = arctan(A/(fc + (5fci?)), we find for the channels 
1 and 3: 



Hi,. 



dk r 

2^ 



2kf 



2kF + Vik + SkF)^ + XA a\{k)ai{k) 
v/(fc + SkF^ + X^] al{k)a^{k) . 



The Hamiltonian H2A follows accordingly by replacing 
5kF <-> 2kF- 

We can now distinguish two situations: (i) \5kF\ *C A, 
where the channels j = 1,3 are gapped and therefore 
transport for Ui,2 ^ A is strongly suppressed, and (it) 
|(5fci?| > A, with no gap. Since kp ^ A, transport proper- 
ties involving the channels 2 and 4 will always be perfect. 
For case (ii), the coupling A will thus not significantly 
affect transport, i. e. up to weak perturbations, the quan- 
tized ideal currents are found. However, for case (i), we 
obtain at T = 



h = {e^/2h){Ui + U2) 



(42) 



just as is expected for absolute Coulomb drag. This simple 
calculation shows that sufficiently large 5k p can destroy 
the perfect Coulomb drag. Experimentally, this could be 
achieved by uniformly shifting the chemical potentials of 
the two reservoirs attached to the same wire. 



4 Weakly coupled reservoirs 

How important is the requirement of good (adiabatic) 
contact to the reservoirs assumed in our theory so far? 
Is it still possible to observe the characteristic effects of 
crossed LL transport and Coulomb drag with weakly cou- 
pled reservoirs? To study this point, we assume now that 
reservoirs are coupled to the QWs by tunnel junctions 
with (identical) dimensionless tunnel conductance Tq <C 1. 
In this section, we consider only the limit of zero tem- 
perature and sufficiently large applied voltage such that 
transport proceeds incoherently at all contacts. 

The applied voltage Ui will then split up into three 
voltage drops: First, the part U[ drops at each junction, 
and U'l at the contact between the nanotubes, where Ui = 
2U[ + U'l- The current h injected into the QW from the 
(say, left) reservoir is then 



/i -roGo[/(|aC/{r« 



(43) 



Here the exponent Vg = 1/g applies to the case of an 
end-contacted (spinless single-channel) QW, while Vg = 
[g + 1/5) /2 for a bulk-contacted QW ||]. If the voltage 
drop [/" across the crossing point is significantly higher 
than U[, the correlation effects described in the previous 
section dominate the transport process, and we may use 
the strong-coupling form of the current W, 



GoTb sgn{U[' ± U^){\U^ ± U^\/Tb 
± 



,1/9-1 



(44) 



Since both Eq. ( [43| ) and (Q) must give the same current, 
the condition [/" ^ [/( leads to the estimate 



Ui,U2:>TB[n{aTBY<' 



-1/(1/3+1-1/3) 



(45) 
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Once this condition is satisfied, typical crossed LL effects 
can be observed even for non-adiabatic (weak) coupling to 
the voltage reservoirs. Evidently, Eq. ([45| ) cannot be satis- 
fied for end-contacted QWs. This is because the power-law 
exponent in Eq. (^5|) is always negative (namely —1), and 
with To <C 1 the inequality cannot be satisfied except for 
unreasonably high Ui. However, assuming that we have 
bulk-contacted QWs with g < \/2 — 1 ~ 0.414, the power- 
law exponent changes sign, and then the condition ( ^5|) 
can be fulfilled even for very small applied voltage Ui . 



5 Electron tunneling 

Let us now address the influence of electron tunneling be- 
tween the QWs. For clarity, we shall focus on the case of a 
strict point-like contact {Lc 0) as in Sec. Assuming 
g < 1, tunneling is then irrelevant under the RG and can 
be treated to lowest order in the tunneling matrix element 
t (unless t is very large). In addition, we restrict ourselves 
to small applied voltages, \Ui,2\,Skp <C Tb, at zero tem- 
perature, where the tunneling density of states (TDOS) 
describing electron tunneling into QW a = 1, 2 (at a; = 0) 
in the presence of the electrostatic coupling carries the 
standard power-law suppression factor Q 



p{E) oc 0{E)E^/3-i 



(46) 



Of particular interest will be the effect of Skp = kpi — kp2 
on the transport properties in the presence of tunneling 
(without loss of generality, we consider Skp > 0). 

Tunneling modifies the currents flowing to the left and 
to the right of the coupling point x ~ 0. Under a golden 
rule calculation, these currents can be easily found from 
the rates for tunneling of a p = ± = mover in QW 1 
into a p' mover in wire 2: 



r, 



dE p{E)p{App, - E) . (47) 



Because of the pointlike coupling, there is no momentum 
conservation, and the tunneling electron can change its 
chirality. In Eq. (|47|), 



Skp/g + ppi - pp>2 



(48) 



where ppa is the chemical potential of a p mover in QW 
a — 1,2 (relative to the respective mean chemical poten- 
tial), which is determined by the applied voltages Ui^2 and 
the electrostatic coupling Tb- Simple dimensional scaling 
gives from Eq. (E^) 



pi- 



/\2/s- 

^p Z pp' 



(49) 



For Skp ^ g\Ui,U2\, the tunneling rates become indepen- 
dent of the applied voltages J7i,2- In this limit, tunneling 
rates are large, and therefore tunneling can potentially 
modify the transport behaviors discussed in Sec. ^ How- 
ever, for Skp ^ g\Ui,U2\, the rates are power-law sup- 
pressed, and hence tunneling rates become very small. In 



fact, since the power-law exponent is larger than the one 
in the absence of tunneling, see Eq. (Q), we expect that 
tunneling leads only to subleading corrections for small 
Skp. The role of tunneling in a concrete experiment could 
then be revealed by simply tuning Skp from zero to large 
values. If this leads to a destruction of the crossed LL 
scenario, tunneling matrix elements must be sizeable and 
should be taken into account for large Skp. 



6 Conclusions and applications 

Let us briefly summarize the conclusions and major find- 
ings of our paper. We have studied transport through 
two quantum wires or carbon nanotubes that are coupled 
along a contact length Lc- Assuming that the wires are in 
the Luttinger liquid state with sufflciently strong interac- 
tions, the main coupling mechanism is of electrostatic na- 
ture, and tunneling provides only perturbative corrections 
to the problem. The electrostatic coupling leads to qualita- 
tively new nonequilibrium transport behaviors compared 
to the case of Fermi liquid wires, namely crossed Luttinger 
liquid effects and Coulomb drag. Crossed Luttinger liquid 
effects are characterized by pronounced zero bias anoma- 
lies, dependencies of transport currents on applied cross 
voltages, and resonant behaviors at finite voltages. For 
extended contacts, very pronounced Coulomb drag effects 
arising at low temperatures were found. For short con- 
tacts, the linear transconductance has a maximum at a 
finite temperature and approaches zero as T ^ 0. 

For sufficiently long contacts Coulomb drag can even 
be perfect at low temperatures, with the transconductance 
approaching its largest possible value e^/2h. We have pre- 
sented an exact solution of the full transport problem valid 
for arbitrary contact length Lf, at the special interaction 
strength g — 1/2. In addition, the nonlinear properties 
of Coulomb drag and their relation to crossed Luttinger 
liquid effects were clarified. Our study has assumed adia- 
batic coupling to the voltage reservoirs. However, as out- 
lined in Sec. ^, for sufficiently strong interactions, the 
same characteristic effects are expected for bulk-contacted 
nanotubes, where the nanotubes are in weak (tunneling) 
contact with the leads. In addition, we have shown that 
the tunneling between the nanotubes should cause only 
subleading corrections to the behaviors outlined here for 
small Skp ^ kpi — kp2. However, by varying this quantity, 
the role of tunneling can be easily determined in practice. 
For large Skp, the characteristic correlation effects will be 
washed out by a sufficiently large tunneling matrix ele- 
ment. 

Let us conclude the paper with two possible applica- 
tions. Crossed nanotubes could be used for voltage am- 
plihcation, e.g. by using the circuit shown in Fig. H. The 
circuit equation Ui[l + RG{Ui,U2)] = Uq can be easily 
solved for g = 1/4 and T = 0, with the results depicted 
in Fig. H Clearly, the voltage amplification ratio 81/2/ dUi 
can be tuned to extremely high values. A more complex 
setup built up of three nanotubes coupled to each other in 
a star-like manner is shown in Fig. If the distance Lx 
between the contact points exceeds the length scales v/T 
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Fig. 5. Amplification circuit based on crossed nanotubes. For 
fixed R and Uo it is possible to achieve \dU2/dUi \ S> f. 
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Fig. 6. Voltage amplification ratio for RGo = 1 at T = and 

g = 1/4 (here Tb = 1). 
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Fig. 7. Triangular setup of nanotubes coupled by (local) elec- 
trostatic interactions. 



or v/Ui, transport proceeds in an incoherent way and can 
thus be modelled as a sequence of crossed LLs. The voltage 
drops Ul and [/j arising at the respective contacts obey 
Ui — Ul + , and now effectively play the role of two- 
terminal voltages. This allows for the direct application of 
the results of Sec. ||, and the whole current-voltage char- 
acteristics of such a setup can be obtained. We only point 
to one interesting consequence of this solution, namely the 
possibility of a voltage measurement involving only elec- 



trostatic contact. Setting, say, C/3 = 0, the current Ii is 
seen to vanish once Ui = —KU2, where k is constant and 
assumed to be known from previous measurements. The 
(presumably unknown) voltage U2 can thus be determined 
by tuning Ui to the point where Ii = 0. 

We thank A. O. Gogolin and H. Grabert for discussions, and 
the Deutsche Forschungsgemeinschaft for support under Grant 
No. GR 638/19-1 and the Gerhard-Hess program. 
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